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We investigate the pattern of chiral symmetry restoration of QCD for N c = 3 and TV/ = 2 
at finite temperature (T) beyond the chiral limit. To this end, we employ the instanton-vacuum 
configuration for the flavor SU(2) sector and the Harrington-Shepard caloron for modifying relevant 
instanton parameters as functions of T. The meson loop corrections (MLC), which correspond to 
1/N C corrections, are also taken into account to reproduce appropriate m q dependences of chiral 
order parameters. We compute the chiral condensate as a function of T and/or m q . We observe 
that MLC play an important role to have a correct universality-class behavior of chiral-restoration 
patterns in this framework, depending on m q : Second-order phase transition in the chiral limit, 
m q — and crossover for m q ^ 0. Without MLC, all the restoration patterns are crossover, due 
to simple saddle-point approximations. It turns out that T* = 159 MeV in the chiral limit and 
T c x = (177, 186, 196) MeV for m q = (5, 10, 15) MeV, using the phenomenological choices for the 
instanton parameters at T = 0. 
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I. INTRODUCTION 
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qi The spontaneous breaking of chiral symmetry (SB%S) has been one of the most important and intriguing subjects 

for decades because SB%S reveals complicated structures of the QCD vacuum. Its restoration at finite temperature 
(T) and/or baryon chemical potential (//) can be understood microscopically as the QCD vacuum effect is diminished 

^i and the system changes its nature drastically at a certain [i-T point, as T and/or fi increase. The restoration of chiral 

symmetry can be observed by the changes in the chiral order parameters such as the chiral condensate, the chiral 

Q\ susceptibility, dynamically-generated effective quark mass via SB^S, weak-decay constant and mass of the Nambu- 
Goldstone (NG) boson, etc. In principle, lattice QCD (LQCD) simulation is the most promising method to deal with it, 

>0 although the sign problem has been a huddle at finite /i. To investigate the chiral restoration at finite T and/or fi, there 

have been also many effective approaches such as the QCD sum rules [IHSj, Nambu-Jona-Lasinio (NJL) model 0HS], 
Dyson-Schinger method [H[H], Polyakov-loop augmented NJL model (PNJL) [5HH], gauge-gravity duality model [151 - 
H5] . hidden-local symmetry (HLS) model [161118) . chiral-perturbation theroy (xPT) [l"9l - f21"] . instanton model [22-25], 
functional renormalization-group method |26[ I27j , and so on. 

The patterns of chiral symmetry restoration are intricately linked with the quark mass. At high T, for QCD with 
two-flavor massless quarks, the associated pattern of chiral restoration is belonged to the universal class of 0(4) spin 
model in three dimension, therefore it is second-order transition. However, with the small quark mass, the second- 
order transition is replaced by a smooth crossover. On the other hand, there is growing evidence show that at low T 
the chiral restoration in the chiral limit (m q = 0) is first- order transition. It suggests that there exists a tricritical 
point (TCP) where a line of critical point, so called 0(4) line turns into first-order transition as /i increases and/or 
T decreases. If the quark masses are nonzero, there should be a critical end point (CEP) in the phase diagram of 
QCD at which the first-order transition line terminates and followed by the crossover when T increases or /i decreases. 
The positions of CEP and TCP of the chiral phase transition have been attracting a lot of interest recently [25H54] . 
Moreover, the chiral (scalar) susceptibility, which stands for a response to the explicit breakdown of chiral symmetry 
by nonzero m q and exhibits the pattern of chiral symmetry restoration, is very sensitive to how the QCD vacuum 
behaves with respect to m q . Hence, it is of great importance to study m q dependence for the pattern of chiral 
restoration in a sophisticated manner. 

In this article, we investigate the pattern of chiral symmetry restoration beyond the chiral limit at finite T. For 
this purpose, we employ an effective action (or thermodynamic potential) derived from the instanton-vacuum config- 
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FIG. 1: An one-loop diagram for the meson correlation function. The solid and dash lines indicate the (anti)quark and meson, 
respectively. 

uration 35, 36]. From the previous work for (T, p) = 0, based on an effective chiral action via the instanton-liquid 
model [24], it turned out that the meson loop corrections (MLC), which correspond to the 1/N C corrections to leading- 
order contributions based on mean-field approximations, play a critical role to reproduce appropriate behaviors of the 
chiral order parameters, such as the scalar susceptibility as a function of m q . Moreover, if and only if MLC is applied 
properly, the effective quark mass, which originates from S%SB and relates to the chiral condensate, showed compa- 
rable m q dependences with those obtained by LQCD simulations [24, 37-39 . Because we are interested in computing 
the chiral order parameters in the presence of nonzero m q , it is necessary to include MLC in our calculation. 

In addition to MLC, we also take into account T modifications on the instanton parameters which are average 
(anti)instanton size (p) and inter-(anti)instanton distance (R). In vacuum, their values are estimated phenomenolog- 
ically as p sw 1/3 fm and R ~ 1 fm, comparable with those from LQCD simulations. Following the previous work 25 , 
we exploit the Harrington-Shepard caloron, corresponding to the temporally-periodic semi-classic solution of Yang- 
Mill's action in Euclidean space [ID]. The fermionic Matsubara formula is also used to evaluate the T dependence 
of the chiral order parameters. We compute the chiral condensate as a function of T and/or m q , and observe that 
MLC play an important role to agree with the universality-class behavior of the chiral-restoration patterns in this 
framework: Second-order phase transition in the chiral limit m q = and crossover for m q ^ 0. Without MLC, all the 
restoration patterns are crossover. We also find that T* = 159 MeV in the chiral limit and T* = (177, 186, 196) MeV 
for m q = (5, 10, 15) MeV by using the phenomenological choices for the instanton parameters: R f=a 1 f m and p « 1/3 
fm. These values of T* are sensitive to the diluteness of instanton ensemble ~ 1/R. Within possible deviation of the 
instanton parameters, we can obtain values for T* comparable with those from LQCD simulations. 

This article is organized as follows: In Section II, we briefly introduce an effective action derived from the instanton- 
vacuum configuration, and MLC are added to the saddle-point approximation as the 1/N C corrections. All the 
ingredients discussed in Section II are extended to finite T case. For this purpose, we introduce the Harrington- 
Shepard caloron and fermionic Matsubara formula. In Section IV, we present and discuss our numerical results for 
the chiral order parameters as functions of T as well as m q . The final Section is devoted for summary and conclusion. 

II. EFFECTIVE POTENTIAL WITH MESON-LOOP CORRECTIONS 

In this Section, we first introduce an effective chiral action derived from the instanton-vacuum configuration, for 
the leading order of N c [2 
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where C, N/V , and A denote an irrelevant constant for further discussions, instanton number density (instanton 
packing fraction) being equal to 1/-R 4 « (lfm) -4 , and a Lagrangian multiplier to exponentiate the 2Nf-'t Hooft 
interacton [35], respectively, a indicates the saddle-point value of an isoscalar-scalar meson field, corresponding to 
SB%S, and we ignored other meson contributions here. The value of a will be determined self-consistently by solving 
saddle-point equations in Section IV. Note that this sort of treatment of mesons are equivalent to usual mean-field 
approximations. Tr represents a trace over color, flavor, and Lorentz indices. The momentum- and m q -dependent 
effective quark mass originated from nontrivial quark-(anti)instanton interactions is parameterized as 
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Here, Mo and p stand for effective the quark mass at zero-momentum transfer and average- (anti)instanton size. It is 
worth mentioning that, from the Fourier transform of the quark zero- mode solution, the effective quark mass can be 
written in terms of the modified Bessel functions |36] . However, we will use the parametrized expression as in Eq. pi) 
for numerical convenience. As indicated in Rcf . [36 , as the momentum transfer goes high enough, the parametrization 
of the effective quark mass becomes different from Eq. (|2| and proportional to l//c 6 . We verified that there appears 
only negligible difference from using M k oc 1/fc 6 so that we make use of Eq. fcty for all the momentum-transfer region. 
Note that in principle Mq is a function of m q as well as T implicitly, as will be shown shortly. 

If we consider quantum fluctuations around the saddle-point values of the mesons, we can rewrite the effective 
action with an additional term, which relates to MLC, indicating 1/N C corrections of the effective action [2"4l l3"Tl l3"5] . 
as follows: 
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where T t = (1,75, it, 1x75) relates to the fluctuations from isoscalar-scalar, isoscalar-pseudoscalar, isovector-scalar, 
and isovector-psuedoscalar mesons. By evaluating over the trace, the effective action can be simplified as 
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where we have denoted l/N c leading (LO) and next-leading order (NLO) contributions. Note that we also have used 
a short-handed notation M k = m q + M k . As understood in Eq. Q, the term inside the square bracket in the last 
line indicates a (one-loop) correlation function for the relevant mesons, propagating with momentum q, as shown in 

Figure 111 The coefficient c\ corresponding to each meson is assigned as 



cf = (-1, +1, +3, -3), cf ] = (-1, +1, -1, - 



(5) 



III. EFFECTIVE POTENTIAL AT FINITE TEMPERATURE 

To investigate the physical quantities at finite T, we discuss briefly how to modify the instanton parameters, p and 
R at finite temperature. We will follow our previous work [25] and Refs. [HI HO] to this end. Usually, there are two 
different instanton configurations at finite T, being periodic in Euclidean time, with trivial and nontrivial holonomies. 
They are called the Harrington-Shepard {JD] and Kraan-Baal-Lee-Lu calorons [2JI12], respectively. The nontrivial 
holonomy can be identified as the Polyakov line as an order parameter for the confinement-deconfinement transition 
of QCD. However, since we are not interested in the confinement-deconfinement phase transition in this work, the 
Harrington-Shepard caloron is chosen here for simplicity. 

Here, we would like to explain our strategy to modify the effective action in Eq. Q as a function of T. As in 
Ref. |43) . the quark zero-mode solution can be obtained directly by solving the Dirac equation in the presence of 
the caloron background. By performing a Fourier transform of this zero-mode solution, one is lead to an expression 
for the T-dependent effective quark mass Mq(T). However, in the present work, we choose a simpler way to obtain 
M (T): Since the effective quark mass can be expressed as a function of the instanton parameters, R and p) [36] . 
instead of solving the Dirac equation directly, we modify R and p as functions of T using the Harrington-Shepard 
caloron, resulting in the T-dependent effective quark mass. As will be shown, this modification together with MLC, 
we obtain compatible results with that given in Ref. |43j and other similar approaches. We now write the instanton 
distribution function at finite T with the Harrington-Shepard caloron as follows: 
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FIG. 2: (Color online) Normalized p/po and n/no as a function of T for N c — 3 (left). Mk as a function of T and absolute 
value of the momentum \k\ (right). 



Here the abbreviated notations are given as: 



P = -b\n[A RSPcut }, p = -bln[A RS (R)}, C Nc = 



4.60 e - 1 - 6Sa BsNc 

^(N c -2)\(N c -iyr 



(7) 



w, 



11 



: iv„ 



27 

T 



Nc 



2V?-1 



2 , HJVc - 2A7> JV 

7r , b — -, n= — . 

3 ' V 



(8) 



Note that we defined f3 and /3 at a certain phenomenological cutoff value p cut and (R) ss J?. Actually only /9 is relevant 
in the following discussions and will be fixed self-consistently within the present framework. Ars stands for a scale 
depending on a renormalization scheme, whereas V3 stands for the three-dimensional volume. Using the instanton 
distribution function in Eq. (pi), we can compute the average value of the instanton size, p 2 straightforwardly as 
follows 1331: 
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where v = (b — 4)/2. Substituting Eq. ([9]) into Eq. Q, the distribution function can be evaluated further as: 



d(p,T) =Cp b - 5 exp [-M(T)p 2 ] , M(T) = \a n T 2 
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The instanton-number density n can be computed self-consistently as a function of T, using the following equation: 

n~»M{T) = [Cr(i/)]«, (11) 

where we have replaced NT/V3 — > n, and Tjv) indicates a T function with an argument v. Note that C and can 
be determined easily using Eqs. B and (11), incorporating the vacuum values of the n and p: C s=s 9.81 x 10 -4 
and fj ps 9.19. Using these results we can obtain the average instanton size p as a function of T with Eq. J9JI. 
The T dependences of the normalized p/po and n/nr, are plotted in the left panel of Fig. [21 As shown there, these 
quantities are decreasing with respect to T as expected. However, even beyond T* ss Aqcei ~ 200 MeV, the instanton 
contribution remains finite. 



IV. EFFECTIVE THERMODYNAMIC POTENTIAL WITH MESON-LOOP CORRECTIONS 



To use the effective action in Eq.(4) to study the pattern of chiral symmetry restoration, one needs extend it to a 
finite T case. For this purpose, we make use of the fermionic Matsubara formula as follows: 
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where w n and fc stand for the fermionic Matsubara frequency assigned as (2n + 1)ttT for n £ X and the three 
momentum of the quark. Substituting Eq. (12) into Eq. (El), one can derive an effective thermodynamic potential for 
a unit volume as 
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In deriving above effective thermodynamic potential from Eq. El), we have taken the following approximations: 1) 
For the integral over q^ : we set a cutoff mass A. 2) We also ignore q^ dependence inside the square bracket of the NLO 
contribution for simplicity. In this way, we have an additional parameter A in the present framework. However, as 
shown in Ref. [24 , since the isovector-pseudoscalar meson, i.e. pion dominates the meson fluctuations, it is reasonable 
to set the cutoff A proportional to m v as follows: 
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Note that, in the above equation, we have multiplied a factor po/p to m^ in order to include T dependence of the 
cutoff mass. Moreover, this multiplication factor represents a correct chiral-restoration pattern of m^, i.e. the mass 
of pion, as a NG boson, increases as SB^S restored. In the left panel of Figure [2] this factor increases as T does show 
a reasonable behavior under chiral restoration transition, i.e. m^ gets heavier as chiral symmetry is restored. We 
also have taken into account the third approximation: 3) As mentioned above, the fluctuation from pion dominates 
the NLO contribution. Hence, we ignored all the meson fluctuations except for that from the pion. According to 
these a ppro ximations 1) ~ 3), we obtain the expression for the thermodynamic potential with the MLC contribution 
in Eq. (13). The energies for the quarks are written as E 2 = k 2 + M 2 and E 2 = (k + q) 2 + M 2 . Accordingly, the 



momentum- and T-dependent effective quark mass read: 
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Here, we have ignored the fourth component of the momentum for brevity, hi — > 0, and verified that this simplification 
does not make considerable difference in comparison to the full calculations. Note that p inside M a>b in Eq. (15). 



The value of M is obtained as a function of m q and T by solving the self-consistent saddle-point equation with 
respect to A: <9il c ff/<9A = 0, resulting in 
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where we have introduced a notation £ 2 = fc • (fc + q) + M a M b + (M a M b + M a M b )/2. The relevant functions Fq~2 after 
summing over n are given in Appendix. Note that N/V in Eq. ( 16 ) is a function of T as discussed in the previous 
section. 

We plot the momentum- and T-dependent effective quark mass M a in Eq. (15) for m q — in the right panel of 
Figure [2J As shown in the figure, M a is a decreasing function of |fc| as well as T, according to decreasing instanton 



effect. Thus, the effective quark mass plays a role of a natural UV regulator and signals the chiral restoration with 
respect to T. Similarly, the m q dependence of Mq can be easily computed by solving Eq. ( 16 ) by putting nonzero m q 
in a self-consistent manner. We will see the m q dependence of the effective quark mass in actail in the next Section. 
However, before going further, we want mention that here is one assumption in obtaining Mq as a function of m q : 
p and R are not dependent on m q but only T . We argue that this assumption is reasonable considering that, since 
these instanton parameters, indicating the QCD vacuum properties approximately, do not correspond to external 
sources such as m q , they behave independently with respect to m q . The chiral condensate, (iq*q) is obtained from 
the thermodynamic potential as follows: 
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Here, E% = k 2 + 
can refer to Ref. 



m q and Gq^2 are also given in Appendix 



For more details of derivation of Eqs. (16) and (171, one 
Since the chiral condensate is one of the most pronounced order parameters for the breakdown 



and restoration of SB^S, we briefly discuss how we study the chiral-restoration pattern using the chiral condensate. 
At certain T, SB%S is restored from the NG phase, following the universality-class pattern of QCD phase transition. 
Thus, at the critical T for the restoration, denoted by T*-, the chiral condensate behaves irregularly. Hence, by seeing 
this irregularity, we can have information on the chiral-restoriation patterns from the order parameters. Another 
way is to observe the behavior of the chiral susceptibility with respect to T [9 . In the present work, following 
Refs. [Till EJ SI] , we sort out T* by differentiating the chiral condensate with respect to T, d(iq^q)/dT. At T*, 
d^q^q) /dT becomes finite maximum for the crossover and infinity for the (first) second-order phase transitions. 



V. NUMERICAL RESULTS 



In this section, we present the numerical results for various quantities, such as the effective quark mass and the chiral 
condensate as functions of m q and T. To identify the effects from the MLC contribution, we will show those quantities 
with and without MLC, separately. In Figure |3l we depict the effective quark mass with zero-momentum transfer 



Mq, which was determined by solving the saddle-point equation in Eq. (16), as a function of m q for T = (0 ~ 150) 



MeV. Note that the left panel shows it without MLC, while the right one with MLC. As shown in the left panel. Mq 
turns out to be a monotonically decreasing function of m q in the absence of MLC. The value of M without MLC in 
the chiral limit is about 325 MeV and becomes smaller as T increases. This observation is consistent with that given 
in Ref. gS] for the case at T = 0. 

However, the situation changes drastically when MLC is included. At T = 0, the curve of Mq increases until 
m q « 10 MeV then starts to decrease after it. In the chiral limit, the value of M is about 190 MeV. It is about 45% 
less than the value without MLC. We note that this value of Mq is considerably different from that without MLC. 
Actually, this result agrees with with the result of the previous work Refs. [37J EH] • The value of Mq is changed from 
350 MeV to 125 MeV by including MLC in [37] and from 567 MeV to 360 MeV in [38] . The MLC reduces M at 
T = ranged from 37% to 63%. The difference between the results in [37J and [35J is due to the different instanton 
parameters R and p. In Refs. [37J [38], it was also argued that the MLC contribution corresponds to the chiral log 
term which provides considerable contributions to the chiral condensate |46j . Hence, considering that (qq) oc Mq in 
general SBxS pictures, this drastic changes of Mq can be understood in the same way since our present framework 
is essentially equivalent to that of Refs. [371 I3H] ■ We indeed verify that the value of M with MLC is very sensitive 
to the instanton parameters. Furthermore we also verified that the fine tuning of the instanton parameters provides 
only negligible changes in the value of T c x . The shape of the curve for T = is well consistent with the LQCD 
simulation [39] . although the magnitude of the curve is rather different. The difference between our result and the 
LQCD simulation can be understood by that the different renormalization scales between them: 0.6 GeV for ours 
and 1.65 GeV for Ref. |3J5]. Thus, the value of Mq is very sensitive to the choice of R so that one may obtain more 
comparable results with the LQCD data by tuning those instanton parameters. 

As T goes higher, the curve of M changes its shape as well as magnitude as shown in the right panel of Figure [3j 
As expected from the decreasing instanton effects with respect to T, the magnitude of the curves becomes smaller 
being similar to that without MLC. On the contrary, the value of Mq in the vicinity of m q — is more sensitive to 
T in comparison to that without MLC, whereas the curves become relatively flat as m q increases for m q > 10 MeV 
for all T values. In other words, focusing on the region m q < 10 MeV, this behavior indicates that the effect of the 
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3: (Color online) Effective quark mass Mo in Eq. ( |15[ ) as a function of current-quark mass m q for the different values of 
150) MeV without (left) and with (right) the meson-loop corrections. 



explicit chiral-symmetry breaking becomes obvious. This manifestation of the explicit chiral-symmetry breaking at 
small m q can be explained again in terms of the decreasing instanton effects, i.e. decreasing nonperturbative QCD 
vacuum effects. Since the strength of the instanton effects prevails over that of the explicit chiral-symmetry breaking 
for lower T, one can see only small difference in Mq with respect to m q , and vice versa for higher T. For the larger 
m q beyond m q sa 20 MeV, this manifestation become obscure since Mq is almost saturated as already mentioned. 
This behavior can be also analytically understood by seeing the functions F\ and F 2 , which contains (E 2 — E?) in the 
denominator as shown in Appendix. In a simple analysis, this term gives 0(m~ 1 ): 



Func(T) 

m 



Func(T) 



El 



(Ak 2 ) ~ 2(AM)(m, + Mo) 



k 2 and the same for AM, and Func(T) is an appropriate function of T. 



(18) 
For Ak 2 > 0, we 



q , enhancing the T dependence, Func(T), 
for T = (0 ~ 250) MeV without (left) and with 



where Afc 2 = k 2 

have AM > 0. Hence, the functions F\ t 2 are enhanced for smaller m 

simultaneously. 

In Figure [4l we show the chiral condensate as a function of m 
(right) MLC. The global behaviors of the curves are very similar to those of Mq. Without MLC, the chiral condensate 
decreases linearly with respect to m q . Again, we observe the manifestation of the explicit chiral-symmetry breaking 
with MLC as seen in the right panel of Figure [4] In addition, for the case with MLC, we see that the slope of the 
curves at m q sa increases then decreases around T = 200 MeV. In other words, at a certain value of T between 150 
MeV and 250 MeV, there appears an inflection point with respect to T. Considering that this slope corresponds to 
chiral susceptibility [23], the inflection point indicates the chiral phase transition at T*. We will see this in detail in 
what follows. 

The numerical results for the chiral condensate are drawn in Figure [5] as a function of T for m q = (0, 5, 10) MeV. 
Since the empirical values for u and d quarks lie below about 10 MeV, m u — (2.55l°]o5) MeV and m d = (5. OCt^H) 
MeV for instance [17], we confined ourselves to small m q values. In the left panel of Figure [5) the chiral condensate 
without MLC shows smoothly decreasing curves, and the magnitude decreases for larger m q as expected from the left 
panel of Figure [4] These smoothly decreasing curves indicate that the chiral restoration without MLC is crossover 
with different m q values. This observation is far different from the universality consideration of the chiral phase 



transition 5, 48 , which is second order for 



and crossover for ra q ^ 0. This flaw is cured by inclusion of MLC 



as shown in the right panel of the figure. The solid line in the right panel of Figure [5] indicates the chiral condensate 



for m n 



0, i.e. in the chiral limit, whereas the dot and dash lines those for m q = 5 MeV and 10 MeV, respectively. 



Obviously, in the chiral limit, we observe the second-oder chiral phase transition, and it becomes crossover for m q =/= 0, 
just satisfying the universality class of the chiral phase transition [JS]. Moreover, the differences between the curves 
are negligible for the region below T < 100 MeV in contrast to those without MLC. This behavior is consistent with 
that shown in the right panel of Figure HI The chiral condensate is insensitive to m q below T w 100 MeV. We note 
that these results are consistent in principle with those given in Ref. [S], in which NJL model was exploited beyond 
the mean-field approximation. However, here we would like to give some comments concerning this universality-class 
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FIG. 4: (Color online) Chiral condensate (iq^q) in Eq. |l7| 

temperature, T — (0 ~ 250) MeV without (left) and with (right) the meson-loop corrections 



as a function of current-quark mass m, for the different values of 



pattern of the chiral restoration via the present framework with the MLC: 1) Although the MLC contribution is 
certainly very important and plays a dominant role among the 1/N C corrections, there are still other 1/N C corrections 
which are not included in our calculation. Consequently one cannot exclude the possibility that the result may be 
different with the full consideration including the complete 1/N C corrections. 2) We modified the instanton parameters 
using the the Harrington-Shepard caloron, which ignores the dynamical quark contributions. The outcome would be 
different if the different types of caloron solutions such as the KvBLL one are used, (see, Refs [5UH52] for recent 
developments based on the calorons with the quark determinant.). 3) We replace the fourth component of the pion 
momentum q by the cutoff mass A oc m^ as in Eq. ( 14 ) . Although we verified that the change of the value of the A does 



not cause qualitatively significant changes in the relevant physical quantities, the T dependence of the A (A oc po/p) 
would be too simple to be real. Even it satisfies the T dependence of the m n : the m^ increases as the SB%S is 
partially restored. Naturally a more realistic T dependence of the A would produce different outcome. Admittedly 
our present results may be altered to a certain extend due to these ingredients mentioned above. Nevertheless, since 
either the full 1/N C calculation or the similar calculation using the KvBLL caloron is formidable task which requires 
much greater effort and time to achieve. On the other hand, our relative simple model calculation did grasp the effect 
of the MLC contribution which is believed to be dominant one , we consider that our current result is still helpful and 
very instructive for the further development of understanding the phase structure of QCD. We would like to leave the 
more sophisticated approach for the future works. 

The critical T for the chiral phase transition, T* can be obtained by differentiating the chiral condensate with 
respect to T as mentioned in the previous Section, resulting in T x — 159 MeV in the chiral limit, as long as the 
phenomenological choices of the instanton parameters R ~ 1 fm and p ~ 1/3 fm at T = are adopted. As for 
rn q = (5, 10, 15) MeV, we obtain T c x = (177, 186, 196) MeV, respectively. According to this observation, T x increases 
as m q does. It is worth mentioning that, from LQCD simulations in the chiral limit, it turned out that T x « 180 MeV 
for Nf = 2 using the clover- improved Willson fermions [53j . Also, using the renormalization-group (RG) improved 
action, it was found that T c x sa 171 MeV [51]. In Ref. [32], employing the Nf = 2 staggered fermions, the critical T was 
estimated as T x = (165 ~ 181) MeV. From effective-model approaches, Ref. [7], using the Schwinger-Dyson model, 
obtained T x in a wide ranges depending on model parameters: T x = (82 ~ 169) MeV. Similarly, the NJL model 
analysis showed T x w 125 MeV with the 1/N C next-to-leading order (NLO) computation [5J. In their work, it was also 
shown that T x appears at higher T without the NLO contribution, i.e. in the mean-field approximation. However, 
even in the case without the NLO contribution, the pattern of chiral symmetry restoration presents second-order 
phase transition, being different from ours. 

Comparing with the LQCD values, our current value (T x = 159 MeV) presents only about 5 ~ 20% deviation. 
Hence, we can conclude that, qualitatively, the assumptions discussed in Sections IV and V can be justified. We also 
note that T x can be much closer to the LQCD values by tuning the instanton parameters within possible changes of R 
and p as in Ref. [37] 155] 155] . We obtain T x = (168 - 177) MeV by choosing R = (0.98 - 0.96) fm. It is also observed 
that T x is very sensitive to R. In other words, it is sensitive to the diluteness of the instanton ensemble. Technically, 
this tendency can be understood by that R controls the instanton-number density N/V ss 1/R , which determines 
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FIG. 5: (Color online) Chiral condensate (iq q) in Eq. (171 as a function of T for the different values of m — (0 ~ 10) MeV 
without (left) and with (right) the meson-loop corrections. 
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FIG. 6: (Color online) Chiral condensate {iq ] q) 1/3 [MeV] in Eq. (17 1 as a function of T and 



all the strength of relevant quantities in the present framework as in Eq. (16 1. Theoretically, it is interesting to see 
that T* is inversely proportional to R as shown above. We can explain this as follows: If the instanton ensemble is 
relatively denser, corresponding to smaller R, the vacuum effect tends to remain effective up to higher T. In contrast, 
in the dilute instanton ensemble, the chiral NG phase which is governed by the vacuum effect will be terminated more 
easily at lower T. In this way, we can understand the higher (lower) T* for the dense (dilute) instanton ensemble, 
respectively. Numerical results are summarized in Table [T] 

Finally, in Figure [6| we plot the chiral condensate as a function of T and m q . Although we have been interested 



only in the SU(2) light-flavor sector, we show the numerical result up to m q « 200 MeV to see overall behaviors 
for the m q dependence of the chiral condensate. It is clear that the second-order chiral restoration at m q — turns 
into that of crossover gradually as m q increases. In addition, the magnitude of the chiral condensate decreases as m q 
becomes more massive. This is a consequence of general behaviors of the QCD vacuum contribution which decreases 
as m q increases [57]. The decreasing instanton effects make the chiral condensate smaller as T gets higher as already 
shown in Figure. [5] 
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VI. SUMMARY AND CONCLUSION 

We have studied the pattern of chiral restoration for the SU(2) light-flavor sector at finite T beyond the chiral limit. 
We used the instanton-vacuum configuration in the grand-canonical instanton ensemble. The relevant instanton 
parameters such as R and p were modified at finite T by using the Harrington-Shepard caloron. As have been 
noticed in the previous works, MLC, which correspond to the 1/N C corrections, are critical ingredients to incorporate 
appropriate m q dependences of the physical observables. Hence, we employed MLC in the present work to investigate 
the pattern of chiral symmetry restoration as a function of m q . We computed the effective quark mass and the 
chiral condensate as functions of T and/or m q . The critical T for the chiral-phase transition, T* was obtained by 
differentiating the chiral condensate with respect to T. Our main results are listed below: 

• As a general consequence in the present framework, the instanton contribution is weakened as T increases 
because the instanton ensemble becomes dilute. It results in the decreasing magnitudes of the chiral order 
parameters. As a result, the magnitude of the effective quark mass as well as the chiral condensate, as the order 
parameters for the chiral-phase transition, get diminished with respect to T. Inclusion of MLC does not change 
this trend. 

• In contrast, their m q dependences are largely dependent on whether MLC is taken into account or not. Without 
MLC, the chiral order parameters are monotonically decreasing functions with respect to m q . If and only if MLC 
is taken into account, the computed effective quark mass can be comparable with that from LQCD simulations. 

• MLC is indeed responsible for appropriate chiral-restoration patterns, which follow the universality-class restora- 
tion patterns depending on m q : Second-order and crossover phase transitions for m q = and m q > 0, respec- 
tively, for Nf — 2. Without MLC the chiral-restoration patterns are all crossover, regardless of the value of 
m q . 

• Our results show that T* = (159, 177, 186, 196) MeV for m q = (0, 5, 10, 15) Mey for the SU(2) light-flavor sector, 
when we employ the phenomenological choices for the instanton parameters, R s=s 1 fm and p ~ 1/3 fm at T=0. 
We observe that T* is sensitive to R, and the LQCD compatible values of T* can be obtained by choosing the 
instanton parameters, which are deviated by a few percent from the phenomenological ones. 

Overall, we conclude that MLC is crucial to produce correct patterns of chiral restoration for Nf = 2 in the present 
framework. This observation is very important when one tries to apply this framework to study the QCD phase 
diagram, in particular, the positions of CEP and TCP. The extension of the present results to the nonzero p and/or 
Nf = 2 + 1 cases are urgent tasks to explore realistic understanding of QCD in medium. The related works are 
underway. 
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Appendix 



The relevant functions in Eqs. (16) and (17) are defined as follows 






-oc 



wl + El 2E a \2T 



Phase transition 



MeV 
Second order 



159 MeV 



5 MeV 
Crossover 



177 MeV 



10 MeV I 15 MeV 
Crossover Crossover 



186 MeV 196 MeV 



TABLE I: T* and phase-transition pattern for different values of current-quark masses. 
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Fi = T 



F 2 = T 
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g 2 = t j: 
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where the definitions of M ab , M a ^ b and E ab are given in the text. 
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